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Abstract 

We use a novel method to calculate analytically the QCD instanton prefactor due to a quark 
field carrying a small mass parameter m. In the SU(2) instanton background of size p, the spinor 
effective action (in the minimal subtraction scheme), which gives rise to the prefactor exp(— F^), 
is shown to have the small-mp behavior 

F^ = -ln(m//i) - ln(^p)/3 - 2a(l/2) - (m/9)2{ln(mp/2) +7 + 1/2} 

-2{mp)^{- ln2(mp)/4 + ln(m/9)(l/2 - 7 + ln2)/2 + C} + 0{{mpf), 

where 7 = 0.577216 . . ., a(l/2) = 0.145873 . . ., and our numerically evaluated value for the constant 
C is C = —0.382727 .... A good agreement between this form and the numerically exact calculation 
is found if (mp) ^ 0.8. 
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I. INTRODUCTION 



Fermion determinants in the background of certain nontrivial gauge fields play some 
important roles in connection with various quantum field-theoretic studies, but they are 
notoriously difficult to calculate analytically. At the heart of the problem lies the fact that 
this requires the knowledge on the spectral shifts for given partial differential operators with 
the introduction of background fields. But, for partial differential operators which are com- 
pletely separable (as would be the case with generic radially-symmetric background fields), 
it has become possible to evaluate the corresponding functional determinants numerically 
with high precision. See for instance Refs. [l, 2, 3] for a particularly efficient method devel- 
oped recently by us (with G. Dunne). Our method, which is really a combination of analytic 
and numerical schemes, was first used to compute the QCD instanton determinant for an 

n 

arbitrary value of quark mass 111. Then, applications to other cases involving more general 

nnn 

classes of radial background fields were also made [2|, |3|, |J] . 

The single- instant on determinant due to a quark field of mass m is essentially a function of 
mp where p is the instanton size (as the dependence on the renormalization mass scale p can 
easily be separated). Aside from finding the precise numerical value for the determinant for 
given values of m and p, it will then certainly be desirable, for practical physics applications 
for instance, to have an explicit (analytic) function form for the quantity in some limiting 
cases at least. For large mp a systematic large mass expansion [5] (accounted for most simply 
by using the heat-kernel expansion) provides such analytic formula, and for mp ^1.3 this 
formula was shown to be a good representation of our numerically determined results. More 
interesting is the small mp-limit, as instanton effects are believed to be important for light 
quarks especially, 't Hooft, in a classic paper j^, first obtained the analytic expression 
appropriate to the mp limit of the determinant and the result of being essentially 
exact for all mp values, is of course consistent with his expression. But, as we increase the 
value of mp from zero, a small discrepancy was observed between the numerical result of j| 

riQ 

and the formula based on the small-mass expansion of the determinant p, 13] given in earlier 
iteratures. As the latter expansion (utilizing the known exact massless particle propagator 
sl in the instanton background) required rather delicate manipulation for certain integrals 
appearing due to the infrared problem, we had naturally our suspiction on this small-mass 
expansion. We thus decided to develop an entirely new approach for the small-mass behavior 
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of the instanton determinant, which is the present work. Here we not only go beyond the 
previous finding but also introduce a new method which should have useful application in 
dealing with other related problems. Our new small-mass expansion formula is now fully 
consistent with the numerical result of 

In this paper we will discuss small-mass behaviors of the instanton determinant within 
the general strategy of Ref. but now utilizing the (analytically developed) small- mass 
expansion for the very part reserved for numerical evaluation in Ref. i.e., for low-partial- 
wave radial determinants. Note that these radial determinants are readily determined once 
one has the global (i.e., valid over the entire radial range) solutions to the Gel'fand-Yaglom 
initial- value problems of Sturm-Liouville-type differential equations [9]. So the nontrivial 
part of our analysis concerns the development of globally- valid small-mass perturbation the- 
ory for the solution of the Gel'fand-Yaglom initial value problem. In view of the presence of 
two characteristic length scales, i.e., p and ^, we here find such globally- valid approximation 
by combining one perturbative solution valid in the radial range < r ^ R (here, R can 
be any value satisfying the condition p <ti R <^ ^) with another perturbative development 
valid in the asymptotic range r ^ R, through an appropriate matching procedure. (A similar 
approach has been employed to discuss the 't Hooft-Polyakov magnetic monopole solution 
for small, but nonzero, Higgs mass in Ref jlOl]). 

This paper is organized as follows. In Sec. [TTl we give a brief review on the calculation 
scheme used in [l, 2, 3| and also collect useful results for our later discussions. Then, in Sec. 

we present a detailed discussion on how the desired small-mass perturbative solution 
to the Gel'fand-Yaglom initial value problem can be secured by combining perturbative 
developments derived in distinct radial regions (through matching them in the region of 
common validity). These results are then used to determine small-mass behaviors of the 
QCD instanton determinant explicitly up to 0{{mpy) in Sec. [IVl In Sec. |V]we discuss our 
findings against the numerical result of Ref. jl| and conclude with some relevant remarks. 
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II. REVIEW OF SOME RELEVANT RESULTS 



instanton field 



Because of a fiidden supersymmetry [6|, tfie spinor Dirac operator in an SU(2) single 



12| 



M^) = ^iP^, (r - v^) (2.1) 

has the same spectrum (except for zero modes and an overall multiplicity factor of 4) as the 
corresponding scalar Klein-Gordon operator. Therefore the one-loop effective action, i.e., 
minus the logarithm of the functional determinant for a Dirac spinor field of mass m (and 
isospin i), r^(y4; m), is directly related to the corresponding scalar effective action T^{A; m) 
(for a complex scalar of mass m and isospin |) by |5| 

r^(A;m) = -ln 2r^(A;m), (2.2) 

/i 

where /i is the renormalization scale. The renormalization prescription for the effective 
actions here is that of minimal subtraction. Thus, for any given mass value m, it suffices to 
consider the scalar effective action to learn about the corresponding fermion effective action. 
Also useful is the fact that T^[A] m) can be related to another quantity V^{mp), which is a 
function of mp only, by 

r^(A;m) = ^ln(/ip) + P(mp). (2.3) 
o 

(The factor | in (12.31) is determined by the one- loop /3-function) . Then, concentrating on 
the mp dependence of V^{mp), we can set the instanton size to unity, i.e., p = I- 

To determine the effective action T^{A; m) in the single instanton background ( 12. ip . one 
should be able to evaluate individual partial wave contributions 

^ , , / det(Hj + m2) \ 

Here, J = (/ = 0, |, 1, |, ■ ■ ■; j = /±|) represent quantum numbers needed to 

designate given partial waves, and 'Hi^ij), Tif'^'^ are radial differential operators given by 

4/(/ + l) + f , 4(j-/)(j+/ + l) 3 
^('.i) = -7^ + -2 + 



Qy,2 J.2 + 1 (r^ + 1^ 

92 



5^2+^efr(r), (2.5) 

^free _ , m + 1) + I 6) 



respectively. (Our notations here slightly differ from those of Ref. Ti-^ij) and Hf^^ of 
the present paper do not represent the covariant Laplacian —D^ in the given partial wave 
(as in Ref. jl|), but the corresponding differential operators after extracting a measure 
factor r'^/^, i.e., coincide with the operators 'H{ij) and Tif^^ of Ref. [2I). Such knowledge on 
individual partial wave contributions may then be used to determine r^{m) with the help 
of the formula (first derived in Ref. [l|]) 

L 

)(m)} + r;>L(m) 



lim 

L— >oo 



Y: {21 + 1)(2/ + 2) {r(,,^.^,+ i)(m) + 



1=0,1- 



(2.7) 



where the partial wave contribution to the effective action, j)(m), is defined by (12. 4p . The 
large partial wave contribution (with the renormalization counter-term included), ri-^L{m), 
was calculated to have the following explicit form: 



1 m 
6 ^ ^ 



2^ , ^ 127 

nL H 

72 



In 2 



In — + 1 

4 



(2.8) 



To achieve faster convergence for L — 00 in (12.71) . one can also include higher WKB terms 
in r;>L(m) as discussed in {2, 3]; but, we have no immediate use for such improvement terms 
in this work. 

To evaluate the partial wave contribution r(/j)(m) for / < L, the method due to Gel'fand 



and Yaglom [9| can be used. That is, given two radial differential operators Aii = 'H{i^j) +m 
and M.2 = 'Hf'^'^ + m? on the interval r G [0, 00), the ratio of their determinants is given by 



detA^i 

— TT- = 1™ 
detAl2 ^^00 



khee I 



R) 



where 0(r), 0^'''^'^(r) satisfy the initial value problems: 

+ m^) 0(r) = 0, ^{r) - 

and 



2/-1-- 
r 2 as r 



(2.9) 



(2.10) 



free ^ ^2 



I free/ 



0, 



I free f 



2/4-- 

r 2 as r 



0. 



(2.11) 



With nonvanishing mass m the analytic solution to (12.101) is not known, and hence in Ref. 
[l| an extensive numerical analysis was necessary to find 0(R) when R is large. With m = 0, 
however, the explicit solution to the given initial value problem for general l,j is found in 
terms of the hypergeometric function 



[0(r)] 



m=0 



Vr^ + 1 



(2.12) 



and consequently we have 



«(0 - l*o(r)]„,,.„., = ^;,^,)^ (2-13) 



r,(r) ^ \Ur%^U,v, = (2.14) 
On the other hand, the solution to (12.111) for any mass value m is 

^free(^)^L + i: _ ^^/^.^.(^r), (2.15) 

V2 Vm/ 
and hence, with m = 0, the expression 

0fr-(r) =r2'+f. (2.16) 

If one uses the results provided by (I2.13p . (12.141) and (12.161) with our formulas (12.91) and 

(12. 7p . the massless limit form of 't Hooft for the instanton determinant can be obtained. 



III. MASS PERTURBATION WITH THE GEL'FAND-YAGLOM INITIAL 
VALUE PROBLEM 

In this section we will find the asymptotic behavior of the solution to (I2.1UI) when m is 
small. To this end, the small-m power series development for (^{r) may be considered, with 
the exact m = solution given in (I2.13P or (12.140 as the zeroth order solution. A naive 
construction of such series works fine as long as r is taken to be finite, i.e., is applicable in the 
range < r ^ i? where p(= 1) ^ i? ^ ^. We will refer to this series solution as the small-r 
solution below. The small-r solution breaks down if r becomes comparable to ^, and hence 
it cannot be used to find the asymptotic form 0(R) (needed to apply (12.91) ). To overcome 
this, we will introduce an alternative perturbation series in m, which we will call as the 
large-r solution. The range of application for the latter series is for any r > i? with i? ^ 1. 
We will then connect this large-r solution to our small-r solution in the intermediate region 
(represented by r = i?) i.e., match the small-r and large-r solutions to obtain the globally 
valid solution to our initial value problem. The asymptotic form 0(R) can be determined 
from this global solution. Finding a global solution in this way of matching is known as the 
boundary layer method. For generic mathematical discussions about this method, see |l^. 
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First we describe a general method to construct the small-m perturbative solution to the 
initial value problem associated with an differential equation 

[-^ + yir?jf{r) = 0, (3.1) 

where V{r) can have a series form 

V{r) = Voir) + mVi(r) + mVsW + • • • • (3-2) 

We here assume that the corresponding solution in the m = limit is known explicitly. We 
write the perturbative solution as 

fir) = /o(r) + mViW + m'Mr) + ■■■ , (3.3) 

and then the leading term in the series, foir) = [fir)]m=o, satisfies the differential equation 

(-|^ + ^o(r))/o(r) = 0. (3.4) 

Using the equations of motion (13. ip and (13.41) . we now have the relation 

^{f'{r)foir)-fir)m} = r(r)/o(r) - /(r)/^'(r) 

= {Vir) -Voir)} fir) Ur). (3.5) 

Plugging in the series (13. 3p into (13.50 then yields the recurrence relations 

I- {fk{r)foir) - fkir)m} = ^ Vir)fUr)foir), (A; = 1, 2, ■ ■ ■ )• (3.6) 
dr ^ 

Integrating (13. 6p from a certain 'initial' point r* to r, we find 

f'kir)foir)~Mr)m 

k 

= f Lin foir*) - fkir*)foir*) + T ^ \/,(t)/,_,(t)/o(t). (3.7) 

•^'"* 1=1 

Since the left hand side of (13. 7p is equal to foir)'^-§^.[fkir) / foir)], we can integrate this 
equation once more to obtain the following expression for fkir): 



fkir) = foir) 



\fkir*) 


r du J 


[foir*) J 


r' foiu)'\ 



fkirlfoir*)- fkir*) foir*) 



k 

+ jytY,v,it)fk-mfoit) 



(A; = l,2,---). (3.8) 



These equations may be used successively to find /i(r), f2ir), ■ ■ ■ for r > r* if the function 
foir) and 'initial' values /i(r*), f[ir*), f2ir*), f2ir*), ■ ■ ■ are given. But the uniform conver- 
gence of the above series (13.30 in a specified domain of r is an issue that should be settled 
separately. 



A. The small-r solution 



Now consider our equation (12.101) . i.e., 



r) r^'"'"2 as r ^ 



(3.9) 



with Kff(r) defined in (12. 5p . To find the appropriate series solution 

0(r) = (f)o{r) + m^0i(r) + n/(t)2{r) H 



(3.10) 



by the above method, it is to be noted that we have here Vo(r) = V^Qir), Viir) = 1, and 
Vk{r) = if A; > 2, in the power series notation of (13. 2p . Also, since the initial condition is 
given at the origin, we take r* — >■ and may impose, from (I2.10p . 

Mr*) 



r* ^0 : 



0o(r*) 



0, 



for A; = 1, 2, ■ • • . 



(3.11) 



With this information we can use (13. 8p to find 0i(r), 02 (^); ■ ■ ■ successively. For the first- 
order term, i.e., for 

du 



Mr) = Mr) 



dt Mty 



(3.12) 



/o 0o(m)2 Jo 

with 00 (^) given in (I2.13P or (I2.14p . the integrals can be performed explicitly, to obtain the 
expressions 



<Pi{r) = 0([(r) 



+ 



ln(r^ + 1) 



X 



(/ + 1) 8(/ + l)(2/ + l) (2/ + l)r2 + 2/ + 2 



(2/ + l)r^ 



16(/ + l)(/ + 2) 



2Fi(l,2/ + 4;2/ + 5; 



(3.13) 



01 (r) = 00 (r) 



ln(r^ + 1) 



4(2/ + 1) 8(/ + l)(2/ + l) 
r2{(2/ + 2)r2 + 2/ + l} 



Fi l,2/ + 3;2/ + 4;-r^ 



(3.14) 



8(/ + 1)(2/ + l)(2/ + 3) 

where we have denoted [Mr)] j={i,j =1+1/2) and [Mr)]j=(i+i/2,j=i) by 4>i{r) and (pi{r), respec- 
tively (in accordance with the convention used in (I2.13P and (I2.14p for 0o(r)). Finding the 
next order term, i.e., 



rr ru 

Mr) = Mr) -T^ dtMt)Mt), 



0o(tt)^ Jo 



(3.15) 



is not easy. However, by using the asymptotic series for (f)o{r) and 0i(r) in fl3.15p . we can at 
least obtain the asymptotic (i.e., vahd for sufficiently large-r) series development for 02 (r) 
with some effort. We will describe this calculation below. 

First, from fl2.13p . (12.141) . (13.131) and (I3.14p . the following large-r asymptotic series for 
(po{r) and 0i(r) can readily be deduced: 

'P^^'^)-'^ \2/ + 2 + (4/ + 4)r2 16(/ + l)r4+ /' ^^"^^^ 
<l>oir) = -''"^'(1-^ + 3^ + ■■■), (3.17) 
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2r2 8r4 

2/ + 1 1 / Inr 4/3 + 14/2 + 12/ + 1 



' j 8(/ + 1)(2/ + 3) ^ ^ 1^8(/ + 1)2 ^ 32(/ + l)3(2/ + 3) 
1 / (2/ + 3) Inr / (8/^ + 12/^ - 42/2 - 87/ - 40) \ 
1^16(/ + 1)2(2/ + 1) 64(/ + 1)3(2/ + 1)2(2/ + 3) J ^ 



(3.18) 



.21+?; 



^^'-^2 <^ \- 



1 If Inr 4/3 + 8/2 + / -2 



(/ + 1) r2 \ 4(/ + l)(2/ + l) 16(/ + 1)2(2/ + 1)2 



1/ Inr 12/^ + 8/3 - 33/2 - 35/ - 8\ | 

1^8(/ + 1)(2/ + 1) ^ 64/(/ + 1)2(2/ + 1)2 ^ ^ ^ 

If these series are used with the function inside the w-integral (13.151) . it is not difficult to see 
that its large-r asymptotic behavior (for given J = (/, j = / + |) or J = (/ + |, j = /)) is 
described by the form 

r dt Mt)Mt) = A,r' + (i?2 Inr + A2) r + + + • • • , (3.20) 

where A^s and -B^'s are some constants. By integrating this series, we may then conclude 
that 

du 



Jo (po[u)^ Jo 



^l„4 , {^2,^ B2 2 , ^3,2, 



= -^r^+[^\nr--^ + ^]r' + ^\n'r + A^\nr + Ci + --- , (3.21) 

where Ci is a certain constant. The constant Ci is not given by this consideration which 
relies just on the asymptotic forms of 4>o{r) and 0i(r); but, one can evaluate it based on the 
representation 

du 



Ci = lim 



0o(^^)^ -'o 



dt Mt)Mt) 

(3.22) 



J Ai ^ f B2 B2 A2\ 2 B3 . 2 
- 1 — r +1— Inr ~^ ^ ) ^ ^ r + ^slnr 



as described in the Appendix. By using the large-r form (13.211) with (13. 8p . the large-r 
asymptotic behaviors for 02 (r) can be found. Exphcitly, for sufficiently large-r, we have the 
series 

2/ + 1 



0+(r) 



128(/+ l)(/ + 2)(2/ + 3) 



1 



Inr 



+ 



4/^ + 20/3 + 17/2 - 28/ - 31 



+ 



^32(/ + l)2(2/ + 3) 256(/ + l)3(/ + 2)(2/ + 3)2 J 

(4/^ + 28/3 + 43/2 + 19/ + 1) In r 



In^ r 



\^64(/ + 1)3(2/ + 1) 64(/ + 1)4(2/ + 1)2(2/ + 3) 



„2l- 



11 
' 2 



4/4 + 20/3 + 73/2 _^ ^3g; _^ g5 
1024(/ + 1)3(/ + 2)(2/ + 1)(2/ + 3) 
f 1 



+ 



21 + 1 
2/ + 2 



+ - 



64(/ + l)(2/ + 3) 
/ Inr 
1^"32(/ + 1)2(2/ + 1) 
/ ln2 r 



4/4 + 12/3 _ _ _ 15 \ 
' 128(/ + 1)3(2/ + 1)2(2/ + 3) ; 
(4/4 _ 4/3 _ 29/2 - 23/ - 4) Inr 



V32(/ + 1)2(2/+ 1)2 64/(/ + 1)3(2/ + 1)3 

12/4 + 36/3 - /2 - 58/ - 29 \ ^ 



512(/+ 1)3(2/ + 1)2(2/ + 3) 



(/^O) 



and 



1=0 



11 

f 2 



192 



Inr 5 ' 
^ 128 



3 ln2 r 



32 



+ 



5 Inr 
'64' 



29 
1536 



+ Cr(/ = 0) + 



(3.23) 



(3.24) 



(3.25) 



In (13.231) and (13.241) . C^{1) denote the constants which can be found with the help of the 
relation hke (13:221) . 

The above asymptotic behaviors of (f)f{r) may be compared with the asymptotics of 



')n(r) and 



given in (I3.16p - (l3.19p . Observe that, as r becomes large, 4>k{r) shows a 



power increase in r (with the leading power increasing with k). Therefore, if r becomes 
comparable to ^, the series (13.101) will not be applicable. Thus we may utilize the series 
expansion ( 13.101) for r in the range < r ^ i? where R can be any value satisfying the 
condition 1(= p) < i? < -. 
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B. The large-r solution and matching 



For r > R, we may contemplate on another approximation scheme. First, since r S> 1, 
we can here approximate the potential Ves{r), defined in (I2.5p . by 



4J(J + 1) + | I -4j(j + l) + 4/(/ + l)-3 ^ 

If we consider only the leading term of this potential, our equation (13.91) reduces to the 
massive free equation (12. lip with j instead of I. Then, motivated by the form (12.150 for 
the free solution, we introduce the rescaled independent variable x = mr and consider a 
function (p{x) = (f){r = ^) to recast (13. 9p as 

r 92 



1 + A^Kff (- ) } ^{x) 



m 



0. 



(3.27) 



Now note that, for m ^ 0, we have the expansion 



1 + ^Kfrf- 



m 



1 + 



+m 



4j(j + 1) + 



x^ 



+ m 



-4j(j + l)+4Z(/ + l)-3 



x'' 



4j(j + l)-4/(/ + l) + 6 



+ 



(3.28) 



i.e., the equation (I3.27P is of the very form for which our general perturbation method, 
described in (I3.ip - (l3.8p . can be applied. Therefore, representing the solution to (I3.27P by 
the series 



Lp{x) = ipo{x) + m^ipi{x) + m^(y92(x) + 



(3.29) 



we may use our iterative relations (13. Sp (with the replacements r — > x, r* ^ x*) to determine 
the functions y9i(x), (y92(x), ■ ■ ■ from the given zeroth order function ifQ^x) and appropriate 
initial values given at x = x*. For this it should be understood that 



K)(r) 
V^i(r) 
V2{r) 



Vo(x) 
Vi{x) 
V2{x) 



1 + 



4j(j + l) + ! 



x^ 



-4j(j + l)+4/(/ + l)-3 

x^ ' 
4j(j + l)-4Z(Z + l) + 6 

x6 



(3.30) 
(3.31) 
(3.32) 



What value should we choose for the initial point x*? We here remark that the small-m 
expansion used in (I3.28P is evidently equivalent to a large-x expansion of V(,f[{—)^ and the 
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series development (13.290 would become completely dubious for extremely small x. But, for 
X moderately small so that we still have x ^ m (i.e., for r = R ^ 1), there is no reason 
to suspect the validity of the series expansion (I3.29p . We thus choose our initial point at 
X* = mR (for R satisfying condition 1 ^ i? ^ ^), and determine ip{x) for any x ^ rriR {x 
can be larger than ^) by using the series (13.291) the explicit form of which is found iteratively 
using the relation (13.81) . We can use the function 0(r) = (f{x = mr) to obtain the large-r 
solution to (12.101) . For the zeroth order solution, we here take the (regular) free solution 

^q{x) = (const.) ^^•^j!-"'"^' ^/x hj+iix), (3.33) 

To be general, one may allow foi^) to contain also the singular free solution (involving the 
Bessel function of the second kind, K2j+i{x)). But, in our case, such, more general, solution 
cannot be matched smoothly to our small-r solution, which does not exhibit singularity for 
finite r, and by this reason it is sufficient to assume the form in (13.330 . 

To determine precisely the functions (pi{x), f2{x), ■ ■ ■ by using ( 13.81) . we need to fix the 
overall constant (in (13.330 ) and also specify the initial values ipi{mR), (p[{mR), ip2{mR), 
ip2{'mR), ■ ■ ■. This is where we must demand on the solution if{x) a smooth matching to 
our small-r solution (f){r) at r = R (or, equivalently, at x = mR). That is, we require that, 
for R satisfying the condition 1 ^ i? <^ ^, 

(pQ^mR) + m'^{pi{mR) + m'^{p2{mR) + ■ ■ ■ 

= MR) + m^MR) + m'^MR) + ■■■ , (3.34) 

where (t>o{R), (t>i{R) and 02(-R) are given by the expressions in (I3.16p - (l3.19p . (I3.23l) - (l3.25p 
with r = R. (Note that, our R denoting a generic value satisfying the condition 1 <C -R <^ ^, 
(I3.34P may be viewed as an equation for the 'variable' R). We remark that there is a crucial 
consistency check (to be discussed below) on our matching condition (13.340 . Based on the 
condition (I3.34p . the initial values can then be readily identified: 

^+{mR) = R^'+l 

^-{mR) = 

ifitimR) = 

12 



' 2/ + 1 (2/ + l)m^R^ {21 + l)m4i?^ 

2/ + 2 ^ 8(/+ l)(2/ + 3) 128(/ + l)(/ + 2)(2/ + 3) ^ 



1 + 



m^R^ 
21 + 3 



+ 64(/ + l)(2/ + 3)^"7' 

\nR . 4/3 + 14/2 + 12/ + 1^ 



4(/ + l)m2/?2 + (^8(/ + 1)2 + 32(/ + 1)3(2/ + 3) 



, (3.35) 
(3.36) 



iftimR) = 

1 



im?R'^ 



InR 



+ 



il^ + 201^ + 171^ - 281 - 3V 



32{l + 1)2(2/ + 3) 256(/ + 1)3(/ + 2)(2/ + 3)2 ^ 
1 / InR 4/3 + 8/2 + / -2 



+ • 



(3.37) 



2m2i?2 \^ 4(/ + l)(2/ + l) 16(/ + 1)2(2/ + 1)2^ 
In R 4/4 + 12/3 _ _ _ 15 ' 



32(/ + 1)2(2/ + 1) 128(/ + 1)3(2/ + 1)2(2/ + 3)^ 
2/ + 5 



+ 



(3.38) 



+ 



16(/ + l)m4i?4 

(2/ + 3) InR I (8/4 + 12/3 _ 43/2 _ gy^ _ 40) - 
16(/ + 1)2(2/+ 1) 64(/ + 1)3(2/ + 1)2(2/ + 3) 
In^R 



64(/ + 1)3(2/ + 1) 



+ 



(4/4 + 28/3 + 43/2 + 19/ + 1) In 
64(/ + 1)4(2/ + 1)2(2/ + 3) 



4/4 + 20/3 + 73/2 ^ ^3g; ^ g5 
1024(/ + 1)3(/ + 2)(2/ + 1)(2/ + 3) 



+ 



21 + 1 



cni) 



ip^imR) = R^'+i 



+ 



\nR 



21 + 2 

12/4 ^ g;3 _ 33^2 _ 35/ _ g 

+ 



(3.39) 



+ 



8m4i?4 m2i?2 \^g(/ + i)(2/ + i) 64/(/ + 1)2(2/ + 1)2 

(4/4 _ 4/3 _ 29/2 _ 23; _ 4) ini? 



\n^R 



+ 



^32(/ + 1)2(2/ + 1)2 64/(/ + 1)3(2/ + 1)3 

12/4 ^ 3g/3 _ /2 _ 5g/ _ 29 



512(/+ 1)3(2/ + 1)2(2/ + 3) 



(/^O) 



(3.40) 



and 



(mi?) 



1=0 



3 



+ 



31ni? 



8m4i?4 m'^R^ 
+ 



3\n^ R h\nR 29 

+ 



- 

+ cr(/ = o)U 



y 32 ' 64 1536 ' ' ) ' ]' ^'^'^^^ 

(The initial first-derivative values also follow from these expressions, as R can be regarded 
as a variable). Using (I3.35P and fl3.36p . we can now fix the overall constant in v9o(x) (see 
f l3.33p ) also, to have 

(2/ + l)(2/ + l)! / 2^2H^ 



\f2 \m 
(2/ + 1)! /2\2'+f 



^/x l2l+2ix), 



V2 



m 



2l+l[X). 



(3.42) 
(3.43) 



The consistency check here is that the small-x series of (13.421) and fl3.43p . under the identi- 
fication X = mR, should precisely match the expressions in f l3.35p and f l3.36p . respectively. 
Using (13. 8p for the first order term, we have 



+ 



du 



ipiijnR) 

ipo{mR) ' JmR ^poiu) 



(3.44) 
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where we denoted 



Fi{x) = (p[{mR)(po{mR) — (pi{mR)(pQ{mR) + dt Vi{t)(po{t) 



rriR 



(3.45) 



Here note that the original Gel'fand-Yaglom initial value problem (12.101) does not know 
about our matching point r = R. This then implies that, for ipi{x) (given by (I3.44p ) to be 
related to the solution to the Galfand-Yaglom initial value problem, the right hand side of 
(I3.44P should not depend on the specific value chosen for R, i.e., the i?-dependence should 
disappear from the expression. As it turns out, this nontrivial consistency check is satisfied 
by our form (I3.44p . See below. 

The integral in (13.451) can be performed explicitly in terms of the hypergeometric function. 
Since (pQ{mR) and ipi{mR), given in (I3.35p - (l3.38p . are in the form of the power series of 
mR (with mR taken to be small), it suffices to use the small- argument expansion for the 
hypergeometric function, to see that the /^-dependences originating from the lower endpoint 
of the t- integral are completely canceled by those from ipo{mR) and (fi{mR). Therefore, 
Fi{x) for the two sets of J are given by the i?-independent expressions 

(2/ + l)(2/ + 3) /x\4'+2 „ _ _ 5 



F+(a;) 
Frix) 



4(/ + l)2m3 



X \ 



41 



m 



- j 2i^3 (2/ + 1, 2/ + -; 2/ + 2, 2/ + 3, 4/ + 5; , (3.46) 

(3.47) 



2F3 ( 2Z, 2/ + -; 2/ + 1, 21 + 2,4/ + 3; x^ 



Using these expressions for -Fi(u), the remaining integral in ( 13.440 can be performed, to 
obtain rather complicated expressions involving the Meijer G function llj 



Lp\{x) 



^0 W 



2^ + 3 „3,i 



X 



V 



-1,0, 2/ + 1,-2/ -3 



+ 



2/ + 3 



4(2/ + l)((2/ + 2)!)2 



X 



X (^2/ + 1, 2/ + ^; 2/ + 2, 2/ + 3, 4/ + 5; x' 
2/ I21+2KX) V 2 



^(2/ + 1) - In 



4/2-3 



4(/ + l)(2/ + l) 

Inm , , , ^, , 

y'^{x) + ^J{x 



16(/ + 1)2(2/ +1) 



4(/ + l)(2/ + l) 



(3.48) 



/ 



/ 



X 



1 

'2' 



X 



X 



y -1,0, 2/, -2/ -2 
4' K2i+i{x) „ 3 



4((2/ + l)! 



— 2F3 2/, 21 + -; 2/ + 1, 21 + 2,4/ + 3; x' 

hi+i{x) \ 2 
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z^(2/ + l)-lnf 
4(/ + l)(2/+l) 
In m 



+ 



4(/ + l)(2/ + l)2 



(3.49) 



4(/ + l)(2/ + l) 

where if) is the digamma function. These expressions are i?-independent, the R dependences 
from (fi{mR)/(fQ{mR) being canceled again by those from the lower end of the w-integral in 
( ]3.44p . This nontrivial cancelation supports that our matching procedure works well. One 
may notice that f l3.48p and fl3.49p contain Inm dependences. These Inm terms arise after 
the Ini? terms in (13.370 and (13.381) get canceled by the ln(mi?) factors originating from the 
low end of the w-integral. From the expression (I3.48P and (13.490 . we obtain the following 
results for x — oo: 

-Inf +^(2/ + l) 



4(/ + l)(2/ + l) 
Inf -^(2/ + l) 



+ 



4/^-3 
16(^ + 1)2(2^+1)2 

I 1 



(3.50) 
(3.51) 



4(/ + l)(2/ + l) 4(/ + l)(2/ + l)2 j 

These results will be relevant in Sec. IIVI (To determine the instanton effective action up to 
0{m?), one does not need perturbation terms higher than ipi{x)). 
We now turn to the second order term, which is given by 



+ 



du 



Lp2{mR) 

(fo{mR) ' JmR ^po{u) 



:Fo(u) 



(3.52) 



where we denoted 



F2{x) = ip2{mR)(fQ{mR) — ip2ijnR)v'oijnR) 



+ 



mR 



dt [Vi{t)^o{t)^i{t) + V2{t)^o{tf] . 



(3.53) 



To find the effective action, we have to consider (p2{x) at x ^ oo. It is easy to see that the 
integral in (13.520 converges when x approaches infinity. Thus 



X 



oo 



(p2{mR) 



+ 



(3.54) 



(pQ{mR) JmR v^o(w)2 

Finding the second order term explicitly is not easy. However, using the small-x series for 
ipoix) and <fi{x), we can determine the mR dependences (with mR taken to be small) in 
the quantity inside the square bracket of fl3.54p . This is explained below. (Note that what 
we need for our purpose is just the expression of ip2{x) for x oo). 
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If small-x series expansions for (po{x) and (pi{x) (available from fl3.42p . fl3.43p . (13.481) and 
(I3.49P ) are used in (13.530 . we find that the mR dependence in ^2(0;) again disappear and 
the result can be described in series form 



^ + — LB2 In — + Z^s + - In — + + 



(3.55) 



where D^s and E^^s are some constants. By integrating this series, we may conclude that 

F2(U) 



du 



4x 



' 2 ^ m 4 2 



+ ln2 + inx + C2 + ■ ■ ■ , (3.56) 

2 m 



where is a certain constant. The constant C2 is not given by this consideration which 
relies just on the small-x forms of '^{)[x) and '^\[x)] but, if one wishes to find out its exact 
value, one can try the numerical evaluation based on the representation 



Co 



lim 



du- 
D, 1 



+ 



E2, X E2 

1 In - 

2 m 4 



D2 
2 



Ert ^ r, X 

+ — \n^ — + Ds \nx 
2 m 



(3.57) 



4a;4 

Explicit evaluation of the constant C2 is given in the Appendix. Using (13.561) in (13.541) . under 
the identification x = mR, we find that all R dependences cancel again: so we are assured 
that our matching procedure works fine for the second-order term also. By this procedure 
we are led to the following expressions as x is taken to be very large: 

(3/ + 2) (4/2 + 4/-l)lnm 



In^ m 



32(/ + 1)2(2Z + 1)2 32(/ + 1)3(2/ + 1)3(2/ + 3) 



(P2 {X) 



ifQ {X) 



64(/ + l)(/ + 2)(2/ + 1)3(2/ + 3)2 

j In^m (6/2 + 17/2 + 12/ + 2) In m 

\32(/ + 1)2(2/ + 1)2 ^ 32/(/ + 1)3(2/ + 1)3 

8/3 + 22/2 + 25/ + 12 ^ 

■ C, (/) + C, (/) , 



(3.58) 



256/(/+ 1)3(2/ + l)(2/ + 3) 



(/ ^ 0) (3.59) 



and 











3 In^ m 


Inm 7 


If 2 (x) 


1=0 




1=0 1 


32 


8 768 



+ Cf (/ = 0) + 02(1 = 0)\. (3.60) 



The constants Cf{l) here may be found using the representation in (13.571) . [But, unlike 
Cf{l), C^(/) have some Inm dependent pieces; see (14. 8p - (14.100 below]. In the next section 
these results will be used to find the instanton effective action up to 0{rn}). 
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IV. SMALL MASS BEHAVIORS OF THE INSTANTON DETERMINANT 



In the previous section we derived the asymptotic expressions that the solutions to the 
Gel'fand-Yaglom initial value problems when mass is small. In this section we will use 
those results to determine small-mass behavior of the instanton effective action. Using 
the Gel'fand-Yaglom formula, the sum of the partial wave contributions T(^ij^i^i^{m) and 
r(-;^i j^j)(m) (defined in (12.41) ) can be expressed as 



r(ij=i+i)(m) + r(;+i^^.^;)(m) 

[(Po{x) + m'^(pi{x) + m'^(p2ix) H ][^oix) + ni?Lpi{x) + ■m^(p2{x) H ] 



= + m'rp^ + m'rp^ + Oim"^). (4.1) 
Then, using ^M), fl330D . fl33T|) and fl338|) - fl3:60D . we obtain following results: 

r.-^ln(|±l), (4.2) 

r« = + 3) u 3) 

' 16(/ + 1)2(2/ + 1)2' ^ ' 

(2) _ J ij{2l + 1) + In 2 16/^ - 32/3 _ ^34^2 _ 93/ _ ^2 1 
' ~ j 8(/ + 1)2(2/ + 1)2 ~ 64/(/ + 1)3(2/ + 1)3(2/ + 3) J 
{'^{21 + 1) + In 2)2 ^ (8/2 + 4/ - 3) (V^(2/ + 1) + In 2) 



16(/ + 1)2(2/ + 1)2 64(/ + 1)3(2/ + 1)3 

384/8 ^ 2240/^ + 4560/^ + 2776/^ - 2992/^ - 5702/3 _ 3529/2 _ ^034/ _ 144 
512/(/ + 1)4(/ + 2)(2/ + 1)4(2/ + 3)2 
+C+(/) + Cr(/) + C^l) + C^il), (/ ^ 0) (4.4) 



and 



^(2) _ In^m Inm / 2 11\ 263 (ln2-7)2 3(ln2-7) 

~ 8~ ^ ~8~ V ^ ~ ^ ~ 24/ ~ 4608 16 64 

+C+(/ = 0) + (/ = 0) + C^il = 0) + C^il = 0), (4.5) 

where 7 = 0.577216 ... is Euler's constant. 

For the leading term, i.e., Tf'\ which is appropriate to the zero mass limit, our result 
(14. 2 p was also obtained in [l|. Using the result (14.21) for Ti (and (12. 8p with m = 0) into the 
formula (12.71) yields the zero mass limit of the effective action 



r^(m = 0) = lim 



(2/ + l)(2/ + 2)ln(' 



2 +1\ . nL 127 n2 

—— + 2L2 + 4L + 

2/ + 2/ 6 72 3 
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0.145873... . 



(4.6) 



This is precisely the 't Hooft value [6|. By using the result (14.31) in (12. 7p (together with 
(12. 8p ). we can find the 0{m?) contribution to the effective action. Explicitly, this 0{m?') 
contribution to T^{m) is 



lim m 



^ 4/ + 3 1 

^ 8(/+ 1)(2/ + 1) ^ 2 



In L + In 



m f, m 1 

— On h 7 + - 

2 I 2^2 



(4.7) 



fl, 



While the m Inm term of this expression agrees with the previous result [7|, the numerical 
coefficient of the m? term is different from the value given in 51]. We believe that (14. 7p is 
correct, since we no longer think that the argument used in Appendix B of Ref. which 
is crucial to obtain the old value, is justified. 

If we use also the result (14.41) and (14.51) in (12.71) . we can obtain the O(m^) contribution to 
the effective action. Here we remark that C2{1) contains some Inm dependences, while Ci{l) 
is m-independent. The Inm terms in 6*2(0 originate from our expression for ^i{x), given in 
(I3.48P and (13.490 . Since the Inm terms in ^pi{x) are proportional to ^q{x)^ the integral in 
(I3.57P for the Inm terms is essentially the same as the integral we have in (I3.44p . Thus the 
Inm dependences in 6*2(0 can be found explicitly and so we write 
f 4/3 + 2/2 - 10/ - 7 ?/'(2/ + 1) + ln2 



c+(/) 



32(/ + 1)3(2/ + 1)3(2/ + 3) 16(/ + 1)2(2/ + 1)2 
4/2-3 ^(2/ + l) + ln2 



lnm + C2+(0, (4.8) 



64(/ + 1)3(2/ + 1)3 16(/ + 1)2(2/ + 1) 



lnm + C2"(/), (/t^O) (4.9) 



and 



C2-(/ = 0) = ^ - 7 + ln2) + C^{1 = 0). 

Here, m-independent constants C2(/) have the representation 

F2(u) 



(4.10) 



C2{1) 



lim 



du 



(poiu) 



H — - [ — —mx — ) + — m X + D^lnx 



(4.11) 



4x4 ' x2 V 2 4 2 y ' 2 

where ^(a^) is defined by (I3.53P with the replacement ipi{x) — » <^i(x) (see (I3.48P and (I3.49P 
for (pi{x)), and D^s and E^s are the same constants that appear in (13.551) . Using ()4.4p 
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and (14.51) together with (I4.8p - (l4.10p in (12. 7p . we can obtain the O(m^) contribution to the 
effective action. This Oim'^) contribution to T^im) is 



In^m Inm /I 



I — + ^ - 7 + ln2j + C| , (4.12) 

where the constant C, as derived from (14.41) and (14. 5p . has the value 



{i/j{2l + 1) + In 2)2 ^ (8/2 + 4/ - 3) (^(2/ + 1) + In 2) 



8(/ + l)(2/ + l) 32(/ + 1)2(2/ + 1)2 

263 



+(2/ + 1)(2/ + 2) {c+(/) + C^{1) + C+{1) + C^{1)} 



2304 



^ 384/s + 2240/^ + 4560/6 + 2776/^ - 2992/^ - 5702/^ - 3529/2 - 1034/ - 144 
^ 256/(/ + l)3(/ + 2)(2/ + 1)3(2/ + 3)2 '^^^ 



If we put together our findings, the scalar effective action T^{mp) for small mp (we 
restored the instanton size parameter p here) can be expressed by the formula 

fV.)^a(i).^{..^.-,4} 

+ (1 - , + ,„2) + c} + 0((mrt») . (4.14) 

Then, using this expression with our formulas (12. 3p and (12. 2p . we obtain the small- mass 
behavior of the fermion effective action as given in the abstract. To evaluate the constant 
C, a rather laborious numerical routine is required. (See the Appendix). With some effort, 
we determined this value to be 

C = -0.382727... . (4.15) 



V. DISCUSSION 



In Fig. [T]we present the predictions based on our formula (14.14 ) against the (essentially 
exact) numerical result for the same quantity first provided in Ref. l|. Notice that the curve 
I — the predictions based on our formula (I4.14p considered only up to 0{{mpY) terms — 
matches the numerical result (represented by dots in the figure) closely for mp ^ 0.2, but 
not so for mp larger than 0.2. [The curve II, which is based on the wrong formula suggested 
in Ref. [sl, starts to deviate from the numerical result at much smaller value of mp\. On 
the other hand, the curve III — the predictions based on our formula (14.141) now including 
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r (m p) 
0.15 i_ 




II - 



0.10 



/ 

/ 



0.05 



0.00 



m p 



III 



-0.05 



FIG. 1: The result of the small mass expansion for T (mp), curves I and III, against the exact 
numerical data (dots) from The curve III is based on our formula (|4.14p considered up to 
0{{mp)^), while the curve I contains the expansion terms only up to 0{{mp)'^). The curve II is 
obtained if one uses the wrong formula from [51 . 

up to 0{{mpY) terms — can reproduce quite faithfully numerical results all the way up to 
mp 0.8. In Fig. [2] we have shown this curve 111 together with the predictions due to the 
large-mass expansion (the curve IV), which were obtained earlier in Ref. j^. As one can 
see, entire numerical results (represented by dots) can be reconstructed rather faithfully just 
by connecting the small-mass-expansion-based curve with the large-mass-expansion-based 
curve at the mid-point (close to the value mp = 1) where the two curves meet! 

For more explicit comparison, numerical values of the small mass expansion and exact 
data for several mp values are presented in Table [U Two cases of the small mass expansion 
(one with the terms only up to 0((mp)^) and the other up to 0{{mp)^)) are tabulated. In 
the first case, the values are close to exact values only when mp < 0.1. Note that, in the 
second case, the errors between approximate values and exact values are less than 10^^ when 
mp < 0.1. The errors increase as mp increases but it is still less than 7 x 10~^ even when 
mp = 0.7. 

Finally we note that Fucci and Kirsten recently have developed a similar small mass 
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~s 

r (m p) 



0.15 




m p 



FIG. 2: Plot of our small mass expansion (14.141) (the curve III) and the large mass expansion from 
(the curve IV), together with the exact numerical result (dots) for T^{mp). 



expansion in a different situation when the background space is a generahzed cone 
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TABLE I: Comparision of the small mass expansion and the exact values 



mp 


up to 0{{mpY) 


up to 0{{mp)^) 


Exact value 


Error 


0.01 


0.145662 


0.145662 


0.14566222 


< 10-6 


0.04 


0.143605 


0.143595 


0.14359532 


< 10-6 


0.07 


0.140299 


0.140227 


0.14022767 


< 10-6 


0.1 


0.136281 


0.136039 


0.13603823 


0.000001 


0.2 


0.121366 


0.118924 


0.11890459 


0.000019 


0.3 


0.108978 


0.0999389 


0.09982892 


0.000056 


0.4 


0.103296 


0.0809004 


0.08059098 


0.000308 


0.5 


0.107238 


0.0624694 


0.06193515 


0.000534 


0.6 


0.123057 


0.0446128 


0.04417211 


0.000441 


0.7 


0.152585 


0.0266823 


0.02740940 


-0.000727 



APPENDIX A: NUMERICAL PROCEDURES TO FIND THE CONSTANT C 



In this Appendix we will describe the numerical method used to find the constant C (see 
(I4.13P ). First we have to determine the values of (/) and Cf{l). For Cf{l), we utilize 
fl3.22p and perform the integral given there explicitly for each I value. When 21 is an integer, 
01 (r) (given in (13.131) and (I3.14p ) can actually be expressed in simpler form. For instance, 
for / = 0, 



0+(r) 



01 {r) 



1=0 



1=0 



(r^ + l)Mn(r2 + 1) rl {4r^ + 3r^ - 6) 

1) ^ (2r^ + 3r2 + 2) 



16v^ 

3 

(r^ + 1)2 ln(r^ 



?r2 



(Al) 



(A2) 



and 0f (r) for other / values can be written in similar simple forms as well. With these 
expressions, the integral in (I3.22p can be performed in closed form. For instance, for / = 0, 
we have 

du 



00 (m) 2 Jo 



dt0+(t)0+(t) 



Li2 i-r^) _^ (r^ + 1) (2r^ - 5) \n{r^ + ^) + 



14^4 _ 21^2 ^ 3Q 



64 



192r2 (r2 + 2) 



1152 (r2 + 2) 



(A3) 



du 



00 {u^ Jo 



dt<P^{t)<P^{t) 
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_ 3Li2 (-r2) 2r^ + 16r^ - 21r^ - 6 (-r^ + 2r^ + 4r2 + 1) \n{r'^ + 1) 

~ 64 ^ 384r2 ^ 6^4 ' ^ > 

where Li2(— r^) is the polylogarithm function. Analogous calculation can be done for / = 
|, 1, |, ■ ■ ■ also. Using these expressions with (13.221) . we can determine Ct{l) exphcitly for 
each /-value, to find the results like 

^1^^°^= lS2+^' ^i^Q ^2^ + 3l56' ' ^^^^ 

C-r(O) = , - = ^ , (A6) 

^ ^ ^ 128 128 ^ V2/ 10368 3456 ^ ^ 

To find 02(1), we use (14.111) : but, here, the integral cannot be performed explicitly. 
Therefore, we resort to a numerical method. Since the integral in ( 14. lip diverges when 
X approaches zero, it is better to divide the integration range into two parts. With an 
appropriate choice of the intermediate value u = u*, the small-x series expansion will be a 
good approximation for u < u* (i.e., use the series ( 13.551) for F2{u), including higher order 
terms), and for u > u* the integral can be performed numerically to obtain a finite value. 
By this method, numerical value for C'^(/) can be obtained with some effort. The values for 
/ = and / = |, for instance, are given as 

C'2+(0) = 0.0517859 . . . , C:^ = 0.0145914 . . . , (A7) 
C'2-(0) = -0.0505132. . . , C^i-^ = -0.0294879. . . . (A8) 



,2. 

Using the thus determined values of C^{1) and C^{1) with our formula (14.131) . we can 
determine C numerically. But, since the Z-sum converges rather slowly, it is desirable to 
employ some acceleration technique. Here, to achieve faster convergence for the /-sum, we 
may utilize the expressions of higher-order WKB terms (up to 0{{j)^)) for the quantity 
^i>L{m) (see Refs. for detailed discussions). That is, after performing the /-sum in 

( I4.13P up to / = L, we add a correction term which can be read off from the coefficient of 

for the quantity Ti^L^m), 

10211 



i- ln(4L) + -1 + A (- ln(4L) - M + i- (-1^ ln(4L) + 

L2 I 16 ^ ^ 32/ L3l8 ^ ' 2/ Ln 32 ^ ^1024 



i 



1 f'l /..X 4371 1 f loll / .^ 1685171 , 
With this correction term, the convergence of the /-sum becomes very fast: we can take 
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L = 10 to obtain the accurate numerical result C — 



-0.382727. 
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